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Abstract. New diffraction-limited speckle images of the Red Rectangle in the wavelength range 2.1-3.3 ^im with 
angular resolutions of 44-68 mas ( [Tuthill et al. 2002; ) and previous speckle images at 0.7-2.2 pm (Dsterbart et al 



997 



Men'shchikov et al. 199£) revealed well-resolved bright bipolar outflow lobes and long X-shaped spikes 
originating deep inside the outflow cavities. This set of high-resolution images stimulated us to reanalyze all 
infrared observations of the Red Rectangle using our two-dimensional radiative transfer code. The high-resolution 
images imply a geometrically and optically thick torus-like density distribution with bipolar conical cavities and are 
inconsistent with the flat disk geometry frequently used to visualize bipolar nebulae. The new detailed modeling, 
together with estimates of the interstellar extinction in the direction of the Red Rectangle enabled us to more 
accurately determine one of the key parameters, the distance D » 710 pc with model uncertainties of 70 pc, which 
is twice as far as the commonly used estimate of 330 pc. The central binary is surrounded by a compact, massive 
(M « 1.2 Mq), very dense dusty torus with hydrogen densities reaching hh ~ 2.5 x 10 12 cm" 3 (dust-to-gas mass 
ratio Pd/p ~ 0.01). The model implies that most of the dust mass in the dense torus is in very large particles 
and, on scales of more than an arcsecond, the polar outflow regions are denser than the surrounding medium. The 
bright component of the spectroscopic binary HD 44179 is a post-AGB star with mass M* w 0.57 Mq, luminosity 
Li, m 6000 Z/q, and effective temperature T* ~ 7750 K. Based on the orbital elements of the binary, we identify 
its invisible component with a helium white dwarf with Mwd ~ 0.35 Mq, Lwd ~ 100 Lq, and T W d ~ 6 x 10 4 K. 
The hot white dwarf ionizes the low-density bipolar outflow cavities inside the dense torus, producing a small H II 
region observed at radio wavelengths. We propose an evolutionary scenario for the formation of the Red Rectangle 
nebula, in which the binary initially had 2.3 and 1.9 Mq components at a separation of ~ 130 Rq. The nebula 
was formed in the ejection of a common envelope after Roche lobe overflow by the present post-AGB star. 

Key words, radiative transfer - circumstellar matter - stars: individual: Red Rectangle - stars: mass-loss - stars: 
AGB and post-AGB - infrared: stars 



1. Introduction 

The Red Rectangle is a spectacular bipolar reflection neb- 
ula around an evolved close binary star (also known as 
HD 44179, AFGL915, IRAS 06176-1036). The object has 
been extensively studied for more than two decades (see, 
e.g., references in Waters ct al. 1998| ; [Men'shchikov et "aT 
199S). Recent diffraction-limited speckle images of the ob- 



ject wi:h 62-76 mas re solution in the optical (0.6-0.8 ^m, 



Ostcrbart et al. 1997 j hereafter Paper 1) an d near in- 
frared (0.7-2.2 /im, Men'shchikov et al. 1995 , hereafter 



lar nebula with pronounced X-shaped spikes, implying a 
toroidal distribution of the circumstellar material. No di- 
rect light from the completely obscured central binary 
could be seen. 

New diffraction-limited images of the Red Rectangle 
in the near-IR (2. 1-3.3 /im) with unprecedented angular 
resolutions of 46-68 mas were recently presented by Tuthill 
ct al. (2002|) (hereafter Paper III). The images were recon- 



Paper II) displayed a compact, highly symmetric bipo- 
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structed from the Keck telescope speckle data using the 
bispectrum speckle intcrfcromctry method. The highest- 
resolution images clearly show two bright lobes above and 
below the midplane of an inclined torus or geometrically 
very thick disk. X-shaped spikes along the surfaces of the 
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conical outflow cavities contribute to the intensity dis- 
tribution of the two bright lobes, making them appear 
widened and even double-peaked. A striking feature of the 
Red Rectangle bipolar nebula is its self-similar appearance 
on scales from 80 mas to V and from the red light to at 
least 10 /im, implying that large grains of at least several 
microns in size dominate scattering. 

Only a few two-dimensional (2D) radiative transfer cal- 
culations of the bipolar envelope of the Red Rectangle 
have been published to date. Yusef-Zadeh et al. (1984) 



first simulated the well-known optical images of the neb- 
ula with X-shaped spikes using a Monte-Carlo scattering 
method. Varying the density distribution and scattering 
properties of dust grains, they found that a torus-like 
configuration with a p oc r~ 2 radial density profile and 
with biconical cavities having a full opening angle of 70° 
is able to reproduce the shape of the nebula. These cal- 
culations gave support to the previously suggested idea 
flCohen et al. 1975| ; pVlorris 1981 ; 



Perkins et al. 1981 ) that 



a quasi-spherical envelope with bipolar cavities can repro- 
duce the shape of the Red Rectangle. 



Extending the previous modeling, Lopez et al. (1997 ) 
applied a Monte-Carlo technique in which they were able 
to not only simulate scattering of the stellar radiation at a 
selected wavelength, but also calculate the radiative equi- 
librium temperature and emission of dust. With multi- 
wavelength radiative transfer computations, they aimed 
to constrain the model by comparing it with the ob- 
served spectral energy distribution (SED) and a decon- 



volved 0'.'2 resolution intensity map at 2.2 /im (Cruzalebes 



et al. 1 996 ) . This work represented the first step in the di- 
rection of a more realistic modeling of the Red Rectangle 
capable of explaining a larger set of observational data 
and to reconstruct reliable properties of the object. 



In the previous 2D modeling presented in Paper II , 
we applied a frequency-depend ent ray tracing method 
( Men'shchikov fc Hcnning 1997 ) to construct a detailed 
model of the Red Rectangle consistent with a much larger 
number of observational constraints. For the first time, 
the model reproduced reasonably well the entire SED 
of the Red Rectangle from the ultraviolet to centimeter 
wavelengths and the highest-resolution (76 mas) speckle- 
interferometry images at 0.656 /im, 0.8 /mi, 1.65 /mi, and 
2.2 /im. The extensive modeling allowed us to derive the 
geometry of a compact circumbinary torus-like structure, 
such as the opening angle u> = 70° of the bipolar cavi- 
ties and the inclination angle 8 V — 7° of the symmetry 
axis. For an assumed distance of 330 pc, the model re- 
constructed physical parameters of the object, such as the 
total luminosity L* ss 3000 Lq, the radius R ~ 30 AU 
and the mass M ss 0.25 Mq of the opaque (Ay « 30) 
torus, the density distribution p cx r~ 2 for r < 16 AU 
and p oc r -4 for r > 16 AU, and very large sizes of dust 
particles (a ~ 0.2 cm). 

Although the model describes reasonably well the large 
number of constraints, our new Keck telescope images of 
the Red Rectangle with unprecedented resolutions (44-68 
mas, Paper III) have shown that the model is not consis- 



tent with the longest-wavelength 3.1 /im and 3.3 /im im- 
ages. As we have already demonstrated in Paper III, the 
model predicts a single elongated peak at this wavelength, 
whereas the Keck image clearly displays two lobes divided 
by a dark lane, very similar to the shorter-wavelength im- 
ages. In the model, there is too much emission from the 
hot grains close to the inner boundary of the torus, which 
implies insufficient optical depths. This finding stimulated 
us to recompute the model taking into account the new 
constraints in addition to the old data. 

This paper presents a new detailed study of the Red 
Rectangle based on our previous model ( Paper II ) and on 
the new constraints provided by the Keck telescope images 
( [Paper III] ). In Sect. | we describe our assumptions and 
radiative transfer model of the dusty circumbinary torus. 
In Sect. [| we present the model results and compare them 
with available observational data. In Sect. || we discuss 
the parameters of the Red Rectangle and evolution of its 
close binary. In Sect. |5| we summarize conclusions of this 
work. 



2. Radiative transfer model 



2.1. General formulation 



The detailed 2D radiative transfer model described below 



is similar to our recent model presented in Paper II. This 



new modeling was stimulated by the speckle images of the 
Red Rectangle with the unprecedented, diffraction-limited 
resolutions of 44-68 m as, based on data collected with the 
10m Keck telescope ( Paper III ). The images extend our 
previous series of high-resolution optical and near-IR im- 
ages of the object ( Paper I ; Paper II ) to 3.3 jxm and place 
stronger constraints on models. The new images gave us 
the opportunity to reanalyze the wealth of existing data 
on the Red Rectangle in the widest wavelength range from 
the far-UV to radio wavelengths, using accurate radiative 
transfer calculations. Following Paper III, we denote the 
Keck telescope near-infrared camera filters at 2.2, 3.08, 
and 3.31 /im as Ks, L c , and Lpah, respectively. 

Our aim was to derive from the observations reliable 
information on the structure and physical properties of the 
Red Rectangle and on the evolutionary status of its cen- 
tral binary, HD 44179. We utilized our 2 D radiative trans- 
fer code based on a ray-tracing method ( Men'shchikov fe 
Hcnning 1997 ), which solves the frequency-dependent ra- 
diative transfer problem in axially-symmetric dusty en- 
velopes for a number of dust components and grain sizes. 
A large parameter space was explored in hundreds of mod- 
els by comparing the resulting spectral energy distribution 
and images to the available observational constraints. For 



more discussion of our modeling approach, see Paper II 
Men'shchikov fc Hcnning (2000| ), and Men'shchikov et al 



(2001|) 
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Red Rectangle: Circumbinary torus 
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Fig. 1. Geometry of the circumstellar environment of the close binary HD 44179 (left panel) and three-dimensional 
re presenta t ion of th e mas sive circum binary torus of the Red Rectangle (right panel) as it appears in the near-IR images 
of Paper I , Paper II , and Paper III in projection onto the sky plane. Schematically shown are four regions of the model 
- the innermost dense torus with bipolar cavities (100 AU radius; dark color), the less dense envelope with a p oc r 
density profile (400 AU radius; medium color), the bipolar outflow cavities (lighter color), and the outer extended 
envelope with a steep p oc r -4 density gradient (4xl0 4 AU; the lightest color). The geometry is defined by the opening 
angle of the cavities, ui — it — « 50° (tp w 130°) and the viewing angle, 9 V ss 11°, between the equatorial plane and 
the line of sight. 



2.2. Geometry and assumptions 

The axially-symmetric model torus of the Red Rectangle 
is shown in Fig. 0. The central binary HD 44179 is embed- 
ded in a very large dusty envelope of an axially-symmetric 
geometry defined by the biconical outflow cavities in the 
otherwise spherical distribution of material. For brevity, 
we call such a structure "torus" , although it displays ob- 
vious differences from a mathematician's torus. 

The density distributions p (r) of the massive torus and 
its outer envelope depend only on the radial distance r 
from the center. The density p (r) of the bipolar outflow 
cavities is much lower than p (r) in the central regions 
of the massive torus (r < 800 AU) , whereas in its outer 
parts the outflow is denser than the surrounding envelope. 
Processes of dust condensation or sublimation put the in- 
ner boundary of the dusty torus at i?i « 14 AU, where 
temperatures of dust grains Td rj 1000 K. The location of 
the outer boundary i?2 can be constrained from below by 
a deep coronagraphic image of the Red Rectangle tracing 
the nebula up to 1' from the center (Van Winckel 2001). 
Although the envelope can well be even more extended, 
we have fixed the model boundary at that angular dis- 
tance. In contrast to other mass- losing evolved stars (e.g., 
IRC +10 216, [Men'shchikov et al. 2001]), any more pre- 



cise location of R2 in this object is unimportant: with the 
steep outer density distribution of our model (p oc r~ 4 , 



Sect. 3.3), 99% of the mass is contained in the compact 
dense torus (r < 100 AU). 

As usual, only the transfer of dust radiation was com- 
puted in this work. The gas component, present in the 



model only implicitly, is described by a dust-to-gas ratio 
Pd/p = 0.01 within the dense torus. Instead of assuming 
an arbitrary value, we used constraints from a detailed 
analysis of the binary evolution in the Red Rectangle (see 
Sect. 4.5 for details) to estimate and fix the total mass 
of the circumbinary material at M w 1.2 Mq, therefore 
determining pd/p- Since there are no reliable constraints 
for the outer regions of the envelope, we adopted p^/ p = 
0.004 outside the dense torus, which is close to the usual 
ratio for the mass-losing evolved stars. 

Since the semimajor axis of the spectroscopically ob- 



served component of the close binary (Sect. 4.1) is signif- 
icantly smaller than Ri , we can replace the binary at the 
center of the system with a single star (see also Sects. [2.3| , 
[2.5[ ) for the purpose of the radiative transfer calculations. 
Furthermore, we assume that the plane of the binary's or- 
bit is parallel to the midplane of the circumbinary torus, 
i.e., that the inclination angle of the orbit is i — tt/2 — 6 V 
(cf.Fig.p. 



Extending the previous modeling (Paper II), we used a 
more realistic approach to describe dust grains, the main 
ingredient of the radiative transfer models. Instead of hav- 
ing dust particles of a single (very large) radius, we assume 
that the grains have a wide distribution of sizes. The radia- 
tive transfer code can handle any (reasonable) number of 
dust components, each of them with an arbitrary number 
of grain sizes. In practice, the number of the dust com- 
ponents needs to be kept at a minimum, balanced with 
available observational constraints. Details of the spatial 
distribution of chemical compositions, shapes, sizes, and 
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other parameters of real dust grains in this object are 
poorly known. Large uncertainties of the real properties 
of dust enormously increase the (free) parameter space 
of the models to explore. Since our purpose was to re- 
construct global properties of the Red Rectangle, we had 



to sim plify the dust model (see Sect. 2.6|) and ignore in 
this work the presence of polycyclic aromatic hydrocarbon 



the two distances. Since now we know that LkHa 208 and 
R Mon are of a different spectral type and nature (mas- 
sive pre-main-sequence Herbig Ae/Be stars), this deriva- 
tion does not seem reliable. The same conclusion can be 
made with respect to another estimate of 280 pc ( [Schmidt 



ct al. 1980), which assumed that the binary is visual and 



molecules (PAHs) and of crystalline silicate dust (Waters 
et al. l |998| ) . These are only minor components of the dust 
and our experience shows that they would not alter results 
of the modeling. 

2.3. Effective temperature 

It is impossible to observe the intrinsic energy distribu- 
tion of the central source in the Red Rectangle directly 
since the close binary is obscured by the optically thick 
torus. The radiation of both components of the binary is 
reprocessed by dust grains, which makes our understand- 



ing of what is going on in the system very difficult. Cohen 



et al. (1975 ) classified the radiation source as a B9-A0 III 
star with re 10 4 K. From high-resolution spectroscopy, 



Waelkens et al. (1992 ) favored lower effective tempera- 
tures of T; re 7500 K (or 800 K) and gravity log g = 0.8 
(or 1.4). [Knapp ct al. (1995Q and [Jura et al. (1997j ) mea- 
sured radio emission from the Red Rectangle, concluding 
that it is produced by a small HII region (radius of ~ 17- 
170 AU for our D = 710 pc). |Knapp et al. (1995| ) deduced 
a much higher temperature of at least 2.5 x 10 4 K for the 
ionizing star and Jura et al. (1997) estimated an electron 
density of ~ 10 6 cm -3 in the region. The UV spectrum of 
the Red Rectangle does not show such a hot continuum 
in scattered light, therefore we can safely assume that the 
other star, a cooler supergiant with T* re 8000 K, is re- 
sponsible for the entire SED of the Red Rectangle (see 
also Sect. 

In this work, we have used several models of stellar at- 
mospher es with effecti ve temperatures between 7500 and 
8000 K flKurucz 1993 ). In the final model, described in 
Sect. [| we adopted an atmosphere parameterized by T* = 
7750 K, gravity log g = 1.5, metallicity [M/H] = -3, and 
microturbulent velocity of 2kms~ 1 . The high-resolution 
spectrum of the atmosphere had been smoothed to a res- 
olution of A/ A A ~ 100 and extrapolated to wavelengths 
A > 10 /im. Although our final model's parameters cor- 
respond to a slightly lower value of log g re 1.2, this dif- 
ference does not have any impact on the results of our 
continuum radiative transfer calculations. 



2.4. Distance 

Accurate estimates of distances D are important for reli- 
able derivation of parameters of individual objects. The 
most frequently cited distance to HD 44179 of 330 pc 



(Cohen ct al. 1975) was obtained by assuming that the 
star is a visual binary, that it has the same bolometric 
luminosity as LkHa208 and as R Mon, and averaging 



has identical components of spectral type AO III, whose 
bolometric luminosities are ~ 360 £q_ each. An upper 



limit of 800-1400 pc ( |Waelkens et al. 1992| ) was based 
on the measurements of photospheric abundances, accu- 
rate assessment of the effective temperature and gravity 
(Sect. p. 3] ), and the assumption that circumstellar dust 
produces a 2.7 mag extinction. Although the latter may 
not be reliable, the distance limit see ms reasonable. An 
estimate of 940 pc ( Knapp et al. 1995 ) was based on the 
upper en d of the luminosity ra nge of 2300-7500 Lq sug- 
gested by Waelkens et al. (1992 ) and the bolometric flux of 
the Red Rectangle. Since the luminosity was chosen some- 
what arbitrarily, it implies an uncertainty by a factor of 
1.8 toward shorter distances. Parallax measurements with 
t he Hipparcos satellite gave the value of 2.62 ±2. 37 mas 
(Pcrryman ct al. 1997), which amounts to a distance of 
380 pc with quite large uncertainties making the whole 
range of 200-4000 pc plausible. 

The widely adopted distance of 330 pc is definitely in- 
consistent with the well-established post-AGB status of 
the Red Rectangle. If the value were accurate, the total 
luminosity of the central binary (estimated from the SED 
assuming spherical symmetry of the nebula) would only 
be 920 Lq ( Leinert fc Haas 198S| ) . This luminosity is too 
low for a single post-AGB star (by a factor of ~ 3-10) and 
it seems to be impossible to account for such a difference 
in the framework of the stellar evolution theory for single 
stars. Although binarity does modify stellar evolution to 
some extent, we still believe that the most natural and 
likely cause of this discrepancy is invali d dis tance. 

The results of our modeling (Sect. 3^2) confirm this 
idea, demonstrating that an almost perfect fit to the SED 
of the Red Rectangle can be obtained only if the interstel- 
lar extinction (proportional to the distance) is properly 
taken into account. The interstellar reddening varies with 
direction in the Galaxy, and the variations increase closer 
to the galactic plane. The Red Rectangle is located at a 
galactic latitude of b = — 11?8, in a region without notice- 
able dense interstellar dust clouds. The published photo- 



metric data for the Red Rectangle (Cohen ct al. 1975 



Kilkenny ct al. 1985) show no significant color variations 
with an average B—V = 0.42 mag. The stellar gravity 
and effective temperature correspond to an intrinsic color 
index (B-V) = 0.18 mag flStraizhis 1977| ). As the ex- 
tinction by the dense circumbinary torus is gray, a color 
excess Eb-v = 0.24 mag is due to the interstellar dust. 
This results in an observational estimate of the interstel- 
lar Ay = 0.74 mag, which is close to the value of Ay = 
0.89 mag we obtained in our model fitting of the observed 



SED (Sect. |3.2|). 

From these data one can estimate D toward the Red 
Rectangle using the photometry and spectral classifica- 
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400 600 
Distance D (pc) 

Fig. 2. Observationally determined interstellar visual ex- 
tinction Ay for stars with known distances D in the di- 
rection of the Red Rectangle (within an angular distance 
of ~ 1° from the object). The stars are labeled with their 
respective HD names. The dashed line shows an average 
dependence Ay(D) for our Galaxy. The solid line displays 
the relation used in our modeling, which takes into ac- 
count the well-known fact that within roughly 100 pc 
from the Sun there is almost no extinction. The dotted 
lines indicate the distance and extinction derived in this 
study. Horizontal error bars, drawn only for the Hipparcos 
data, correspond to the uncertainties of the measured par- 
allaxes. The vertical bars are approximate uncertainties 
associated with the color excess estimates. 



tion of nearby stars available in the literature, as well as 
the Hipparcos parallaxes. Our search for the stars was fa- 
cilitated by the SIMBAD and the Two Micron All Sky 
Survey databases. Although not many stars around the 
Red Rectangle have the estimates of spectral types, the 
existing data allow us to put some constraints on the re- 
lationship between Ay and D in this direction (Fig. |J). 
The stars located closer than 350 pc have Ay < 0.3 mag 
demonstrating that a distance of 330 pc (Cohen ct al 
1975|)~]s inconsistent with Ay of 0.74-0.89 mag derived 
in this paper. Data on two B-type stars suggest that Ay 
reaches 1 mag at ~ 800 ± 100 pc, although there is a very 
large gap in the diagram, in the region of our interest. A 
rough interpolation between the few available points seems 
to suggest a distance between 600 and 800 pc, however, 
there are no stars with known spectral types in the region. 
Photometric and spectroscopic observations of V ~ 9-10 
mag A- and B-type stars in the neighborhood of the Red 
Rectangle on the sky would be very useful in obtaining a 
better estimate. 

In this study we adopted an average interstellar extinc- 
tion per hydrogen column number density Ay/N(R) = 



5.3 x 10~ 22 mag cm 2 (Bohlin ct al. 1978) and a hydrogen 
number density njj = 1cm -3 . Assuming also that there is 
no extinction within 120 pc from the Sun, this resulted in 
A v = 1.5 (D/lkpc) - 0.18 mag (Fig. |). Using the value 
of Ay = 0.89 mag derived in this modeling and the depen- 
dence Ay(D) shown in Fig.|^, we can determine a distance 
D w 710 pc, with possible uncertainties of ~ 100 pc due 
to unknown Galactic extinction beyond 350 pc from the 
Sun. 



2.5. Luminosity 

With the above distance of 710 pc, the object's bolo- 
metric luminosity would be 4260 Lq , if the nebula were 
spherically-symmetric. Although the value would better 
agree with the AGB status of the star, this simple estimate 
based on the assumption of spherical geometry is incor- 
rect, because the Red Rectangle is strongly non-spherical. 
Large deviations from spherical geometry in optically 
thick environments usually lead to large errors (by a fac- 
tor of several) , depending on the optical d epths, viewing 
angle, and degree of non-sphericity (see e.g., Men'shchikov 
|fc Henning 1997fc [Men'shchikov et al. 1999| ). O ur present 
modeling shows that the actual total luminosity of the 
binary is L* « 6050 Lq (Sect. [T^). This entire luminos- 
ity seems to be produced by the post-AGB component 
of the binary. Luminosities, masses, and other parame- 
ters of both components, not related directly to the radia- 
tive transfer calculations, are discussed in more detail in 
Sect. |. 



2.6. Dust particles 



As we have shown in the previous modeling ( Paper II ), 
the dominant particles in the dense dusty torus of the Red 
Rectangle are very large. Much smaller amounts of small 
dust grains and PAHs seem to exist somewhere in the neb- 
ula, perhaps in its outer dilute parts and in the outflow 
regions. In the final model we assumed that there are two 
(spatially separate) components of dust grains. Following 



Paper II, the first component is made of very large par- 



ticles of an unspecified chemical composition, with sizes 
of 0.2 cm. This component exists only in the compact, 
extremely dense circumbinary torus, where most of the 
grains are assumed to have coagulated into the large parti- 
cles. The other component is the amorphous carbon grains 
with radii a in the range of 0.005-600 /im and with the 
standard interstellar size distribution dn/da cx a -3 ' 5 . This 
second component exists beyond the boundary of the com- 
pact torus (100 AU < r < R2), where the density drops by 
several orders of magnitude. This particular choice of the 
dust parameters is somewhat arbitrary, as there are no ob- 
servational constraints to better restrict the freedom. We 
should note, however, that it was impossible to make our 
model explain all observations when considerable amounts 
of small dust grains were put in the dense torus. We can 
conclude that a difference between the dust properties of 
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Table 1. Model parameters of the post-AGB star and circumbinary torus of the Red Rectangle. 



Parameter 



Symbol Value Units Uncert. Comment 



Distance 
Total luminosity 
Effective temperature 
Stellar radius 
Inner boundary 
Outer boundary 
Total mass of envelope 
Gaussian density profile 
1st density exponent 
2nd density exponent 
Cavity opening angle 
Viewing angle 



D 


710 


pc 


±10% 




6 050 


In 


±20% 


T* 


7 750 


K 


±10% 


R* 


43 






Ri 


14 


AU 




R-2 


43 000 


AU 




M 


1.2 


Mq 




FWHM 


26 


AU 




(Vl 


-1.5 






a 2 


-4.0 




±20% 


U) 


50 


O 


±10% 


Oy 


11 


o 


±10% 



2.5 



derived (from SED and reddening, Sects, ^t^^) 
derived (from SED and distance, Sec ts 
derived (from observations, Sect. |2.3[ ) 
derived (R t = (Ana)- 1 / 2 ^ 2 T~ 2 ) 
derived (Ri 



3.3) 



71 R*, Sect _ 
assumed (from deep optical images, Sect. 2.2) 
derived (corresponds to Pd/p = 0.01, Sects. 2^. 
model (14-100 AU, Sect.O) 
model (100-430 AU, Sect. |3.3|) 
model (430-43 000 AU, SectTfJ) 
derived (from high-resolution speckle images, Sec t 
derived (from SED and model images, Sects. 3.2, 3 



the compact dense torus, made up of very large parti- (Greaves & Holland 1997; Jura 
cles an i of those of the surrounding low-density envelope, Contreras et al. 1998). To make sure that the observed 



Turner 1998; Sanchez 



made up of normal grains, is a firm result of the modeling. 



3. Results of the modeling 

3.1. Model parameters 

The final parameters of our detailed radiative transfer 
model and their approximate uncertainties are presented 
in Table [|. The model results are described in Sects. 3.2 



3.6 in detail. The density distribution is assumed to be 
a broken power law p oc r tti , the exponents on (i = 1, 
2) being allowed to change in different radial zones. The 
third column lists the values adopted in the final model. 
The fifth column contains our estimates of the uncertain- 
ties associated with the final model parameters. They in- 
dicate very approximately the maximum changes of the 
parameters, such that the fit to the observations can still 
be restored by adjusting some other model parameters. 
The uncertainties should not be interpreted as absolute 
error bars. If some of our general assumptions turn out 
to be insufficiently realistic, this would possibly affect the 
results more than these uncertainties imply. 

It may be useful to recall here that the results of our 
radiative transfer modeling are invariant with respect to 
the distance D, except for several parameters which are 
scaled in the following simple way: 



LocD 2 , A/oc£> 2 , RocD, pcxD" 1 , 



(1) 



where mass and density refer to the properties of dust; for 
conversion to the gas parameters, pd/p must be specified 
(Sects. 2.2, 4.5). These relations make it easy to scale the 
model results to another distance, if necessary. 

3.2. Spectral energy distribution 



In addition to the fluxes used in Paper II, in our new 



modeling we included data collected by the Infrared Space 
Observatory's short-wavelength spectrometer (ISO SWS, 
Waters et al. 1998) as well as new (sub-)millimeter fluxes 



SED of the Red Rectangle is not influenced by obser- 
vations at different epochs, we ignored older broad-band 
fluxes (in the wavelength range covered by the ISO data) 
containing essentially no additional information on the 
spatial distribution of intensity. As our model does not 
contain any species producing uni dentified nebular re d 
emission between 0.55 and 0.75 p,m ( Schmidt et al. 1980 ), 
we included an optical spe ctrum of the innerm ost few arc- 
seconds around HD 44179 ( Reese fc Sitko 199Q that shows 
no extended red emission. 

Figure ^ compares the model SED to the broad- 
band flux distribution and spectrophotometry of the Red 
Rectangle. To emphasize the influence of the bipolar ge- 
ometry, the SED of an equivalent spherical envelope is 
also shown. The latter differs from the torus model only 
by the absence of the outflow cavities (i.e. ip = 180°, ui 
= 0, Fig. |l|). The model of the Red Rectangle fits the 
observed continuum underlying the emission features al- 
most perfectly, from the far UV to centimeter wavelengths. 
Three lower radio fluxes in small beams with a spec- 
tral index of 1.5 are most likely produced by a compact 



H II region ( Jura et al. 1997 ) , whereas three higher fluxes 
in larger beams are mostly due to dust emission. When 
the compact radio emission of the H II region, which has 
not been included in our modeling, is added to the total 
model fluxes, the large-beam fluxes become perfectly fit- 
ted (Fig. |). There is no discernible effect of different aper- 
tures (larger than a few arcseconds) at any wavelengths 
in the model SED, which implies that the dusty torus has 
density distribution highly concentrated toward the center 



(see Sects. 3.3, 3.4). 

The Red Rectangle differs from other dusty objects 
in that the dominant dust particles are very large, most 
likely millimeter-sized (Jura et al. 1997; Paper II). The 



new modeling presented here confirms the idea that the 
dense torus contains most of its circumstellar dust mass in 
large particles having gray extinction. We found no model 
with significant amounts of small grains that would be able 
to reproduce all high- resolution images and the SED. It 
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Fig. 3. Ob served SEP o f the Red Rectangle compared to our model with = 6050 Lq (solid line). The stellar 
continuum ( Kurucz 1993 ), the total model fluxes without interstellar reddening, the model fluxes reddened by the 
interstellar extinction (for D — 710 pc), and SED for an equivalent spherical envelope are plotted. The latter emphasizes 
the major influence of the bipolar outflow cavities on the emerging continuum of an optically thick torus. The observed 
fluxes are shown by different symbols to differentiate between broad-band fluxes, IRAS fluxes, spectrophotometry, and 
ISO SWS data. If available, error bars for the fluxes are drawn, whenever they are larger than the symbols. The model 
assumes that we observe the torus close to edge-on orientation, at 8 V = 11°. Three arrows at the top axis indicate the 
location of the H, K, and L photometric bands. The insert compares the centimeter-wavelength segment of the model 
SED with available radio observations. The dashed line is fitted to the compact emission fluxes from ionized hydrogen 
(Jura et al. 1997). Dash-double-dotted line through the extended fluxes is the sum of the model fluxes (solid line) and 
the compact HII emission. The far-UV model fluxes rising toward shorter wavelengths at A < 0.14 /im and a bump on 
the reddened model curve show the upper limit for the possible contribution of a hot companion with a temperature 



of 6 x 10 4 K (see Sect. 4.1) 



seems to be impossible, however, to fit the observed fluxes 
with only very large grains. 



It was clear also from our previous model ( Paper II ) 
that gray circumstellar extinction alone cannot reproduce 
the SED of the Red Rectangle in the optical and UV. The 
steep peak of the scattered intrinsic stellar energy distri- 
bution of the post-AGB component in the visible sharply 
contrasts the almost horizontal SED in that region. At 
wavelengths shorter than the Balmer jump, the observed 
fluxes drop much deeper and much steeper than simple 
gray scattering of stellar radiation would predict (Fig. ||). 
This puzzle can be solved in a simple and natural way just 
by taking into account interstellar extinction. 



Usually it is difficult to discriminate between the ab- 
sorption by circumstellar and interstellar dust: small cir- 
cumstellar dust grains absorb radiation the same way the 
interstellar grains do. In contrast to the latter, dust in stel- 
lar environments is embedded within and plays an active 
role in modifying the circumstellar radiation field. On the 
other hand, the interstellar dust is just passively absorb- 
ing and scattering photons on their way to the observer. 
Since our modeling has shown that the dust particles in 
the circumbinary torus of the Red Rectangle are predom- 
inantly very large, we are left with only one possibility 
to explain the discrepancy at short wavelengths, namely, 
that the SED is altered by the selective extinction of small 
interstellar dust grains. 
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The distance to the Red Rectangle is large enough 
for the interstellar extinction (Sect. |2.4| ) to affect the ob- 
served SED. For the wavelength dependence of the in- 
terstellar extinction, we used the analytic fits of A\/Ay 
given by Cardclli ct al. (19*89] ) . The extinction law param- 
eter Ry was fixed at the generally accepted value of 3.1 
corresponding to the diffuse interstellar medium. The pro- 
nounced effect of interstellar extinction is visible in Fig. || 
as the increasingly larger difference toward shorter wave- 
lengths (A ^ 1 fim) between the total (unreddened) model 
fluxes and the model fluxes reddened by A\ (Ay ~ 0.89 
mag). After the interstellar reddening has been taken into 
account, the model SED became almost indistinguishable 
from the observed flux distribution. 

In contrast to the transfer of radiation in circumstel- 
lar envelopes, interstellar extinction does depend explic- 
itly on the distance D to the object. The accurate fit of 
the short-wavelength part of the SED enabled us to de- 
rive a distance D w 710 pc to the Red Rectangle, with 
approximate model uncertainties of about 10%, which is 
consistent with both the interstellar reddening data in this 
direction and the object's luminosity. In fact, this larger 
distance implies the stellar luminosity is more than four 
times higher, reducing by the same factor the large dis- 
crepancy between the post-AGB nature of HD 44179 and 
its low apparent luminosity (Sect. ^^)- The actual lumi- 
nosity L+ « 6050 Lq obtained in our modeling is even 
higher than the value of 4260 Lq one would expect for a 
spherical envelope. This is not surprising, given the ob- 
vious non-sphericity of the optically thick circumbinary 



Angular distance 8 (arcsec) 
10- 1 10° 10 1 



torus (Sect. 2.5). 

The above ~ 40% increase of the derived luminos- 
ity is a result of accurate 2D radiative transfer calcula- 
tions. In the optically thick non-spherical configurations 
like the torus of the Red Rectangle, a large fraction of the 
central star's radiation can escape in the optically thin 
polar cavities, thus changing the value of the luminosity 
an observer would derive under the (invalid) assumption 
of spherical symmetry. The observer would underestimate 
the luminosity if the object's orientation is close to edge- 
on, whereas one would significantly overestimate the lu- 
minosity if one sees the object almost pole-on. Corrections 
for converting the observed "spherical" luminosities into 
the real luminosities of the central energy sources may 
reach a factor of several, depending on the distribution of 
opticali deptho, opening angle of the oavitico, and the vie* 



Temperature of large grains 
Temperature of normal grains 
Total (gas+dust) density 
Total density: polar cavities 
Density of large grains 
Density of normal grains 



ing an gle of the object (nee, e.g., Mcn'ohchikov fc Hcnning 



1997 |; pen'shchikov et al. 1999| ). It seems to be difficult 
(if possible at all) to determine such factors once and for 
all; they have to be derived for each individual object in 
a detailed modeling. 



3.3. Densities and temperatures 

The density and dust temperature distributions of our fi- 
nal model of the Red Rectangle are shown in Fig. |i| for 
the smallest grains; the total gas plus dust density is also 
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Fig. 4. Temperatures and densities of dust grains in the 
envelope of the Red Rectangle. The distributions are dis- 
played for the smallest grains of the dust components. 
Also shown are the total (gas+dust) density profiles in 
the midplane and in bipolar outflow cavities along the 
symmetry axis. The total mass of the circumstellar mate- 
rial is 1.2 Mq, corresponding to Pd/p = 0.01 within the 
dense torus. Additional labeling above the lower abscissa 
shows possible expansion times if the outflow velocity were 
constant and equal to 5kms _1 (similar to the observed 
velocities, Jura ct al. 1995| ; Knapp ct al. 2000| ). 



plotted assuming a spatially uniform dust-to-gas mass ra- 
tios pd/p = 0.0 1 within the dense torus and 0.004 outside 
(Sects. 2/2, [4. 5|) . Although pd/p may be a more complex 



function of the position inside the nebula, we assumed this 
simple form in the absence of better constraints. This pa- 
rameter is only a means to approximately convert the dust 
densities into the total density. 

A unique property of the Red Rectangle is that even 
the L c and Lpah images show almost the same intensity 



distribution as the images at shorter wavelengths (Paper I 



Paper II ; Paper III). Moreover, the biconical shape is pre- 



served across a wide range o f wavelengths (0.6-10 /im) and 
spatial scales (0.001-1') (cf . 



Cohen ct al. 1975: Perkins ct 



al. 1981; Waters et al. 1995 ; Paper II; Van Winckel 2001) 



This fact provides a strong observational constraint on the 
density structure and dust properties of the circumbinary 
environment of the Red Rectangle. 

The invariant biconical shape of the nebula at such 
different scales and even at the mid-infrared wavelengths 
suggests that the shape is produced by scattering. The 



10 /im image presented by Waters et al. (1998), in which 
the bipolar outflow can be traced to 10" from the central 
stars rules out thermal emission from dust: in our model, 
grains in the torus and in the cavities have temperatures 
of 32 K and 95 K, respectively, at these distances. At 
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optical as well as at near-IR wavelengths, the biconical 
shape is caused by the radiation scattered by dust grains 
in the bipolar outflow regions. The particles must be very 
large compared to the normal (small) interstellar grains 
since the latter have negligible albedo in the mid-infrared. 
Dominant grains must have sizes of at least several mi- 
crons in order to efficiently scatter the mid-IR radiation. 

The present modeling has shown that the density dis- 
tribution of the toroidal structure of the Red Rectangle 
has two distinct regions: the compact dense torus and 
the extended low-density envelope (Figs. [I], The model 
density distribution displayed in Fig. ^ follows a Gaussian 
profile in the very dense torus, with a full width at half 
maximum of 26 AU and a maximum total density of 
p w 4.2 x 1(T 12 gem" 3 (n H ~ 2.5 x 10 12 cm- 3 ). The 
dense torus contains 99% of the total mass of the cir- 
cumbinary material. The density of the toroidal core drops 
by 5 orders of magnitude at r s=s 100 AU (0715), where it 
changes to a much flatter p cx r~ 1,5 radial profile, which 
in the outermost regions at r > 430 AU (076) steepens to 
a p oc r -4 distribution. Constrained by our new highest- 



resolution speckle images (Paper III, see also Sect. 3.4) 



the radial density profile differs from that of the previous 
model ( Paper 1% by the much steeper, exponential fall-off 
in the outer regions of the torus and by the flatter density 
distribution just outside it. No simpler power-law density 
profile has been found to reproduce the well-defined bright 
bipolar lobes which are almost invariant between the H 
and Lpah bands. 

The bipolar appearance of the Red Rectangle implies 
much lower density in the outflow cones compared to the 
torus over the central subarcsecond region (Fig. ^). The 
inner dense torus is optically thick (r 47 up to A ~ 
1 cm) and our modeling demonstrated that the observed 
appearance of the Red Rectangle requires the outflow cavi- 
ties to be optically thin (t « 0.2). On the other hand, the 
model shows that the outer regions of the circumbinary 
material (at angular distances 6 <; 0'.'5) are optically thin 
in all directions. This implies that the strong biconical 
shape visible in the outer parts of the Red Rectangle is 
produced by dust scattering in the outflow regions which 
are much denser there than the rest of the torus at the 
distances of r > 800 AU from the center. This conclusion 
implicitly assumes that the dust properties in the outflow 
cones and in the rest of the toroidal envelope are the same. 
If dust grains in the torus and in the outflow cones differ 
in such a way that the outflow regions scatter significantly 
more radiation to the observer, there may be no need for 
a density contrast to explain the observations. 

In the absence of reliable constraints, we modeled the 
radial density distribution in the outflow cavities by a con- 



stant density p s 
at r < 2000 AU 



9 x 10~ 19 gcm~ 3 (nn ~ 5 x 10 5 cm~ 3 ) 
(2'.'8) and by a steep p oc r~ 5 profile 
outside (Fig. . The density profile of the outer regions is 
constrained by the long-wavelength fluxes of the observed 
SED, which must be produced largely by the dense in- 
ner torus with very large dust grains. At the distance of 
2000 AU, the density in the outflow cavities is larger by 
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Fig. 5. Mass-loss history of HD 44179 within the torus' 
midplane and in the bipolar cavities, reconstructed from 
the density distribution of the Red Rectangle assuming 
a spherical outflow with constant velocity v — 5 kms . 
Although these assumptions may not precisely describe 
the entire evolution of the envelope, the plot gives an ap- 
proximate overview of the mass-loss history. 

a factor of 14 than that of the toroidal envelope; the den- 
sity contrast is vanishing toward the outer boundary of 
the model. At the inner boundary, the torus is denser by 
a factor of ~ 5 x 10 6 than the outflow cavities. Such a 
density structure can be interpreted as a result of higher 
mass loss and velocity in the outflow along the symmetry 
axis during the last ~4x 10 4 years, as compared to the 
early mass loss that made up the rest of the outer toroidal 
envelope. 

With the assumption of a constant-velocity spheri- 
cally-symmetric mass loss, one can derive the mass-loss 
history of the Red Rectangle from the density profiles 
shown in Fig. ||. The reconstruction presented in Fig. || 
can provide approximate but useful information on the 
global properties and temporal evolution of the mass loss. 
The very high density contrast of the toroidal core of our 
model with respect to the outer envelope (Fig. |j) is re- 
flected in the mass-loss rates that seem to reach 4 x 10~ 2 
AfQyr^ 1 within the dense core, very high values com- 
pared to the usual spherical mass loss from AGB stars. 
However, the dense core is most likely the result of a 
complex, time-dependent hydrodynamical evolution of the 
self-gravitating envelope ejected from the central binary. 
Both expansion and infall are likely to have played a role in 
the accumulation of the circumbinary material and in the 
compression of the compact torus to very high densities. 
For an illustration, note that 90 % of the dense torus' mass 
resides within r = 57 AU. If the mass is redistributed over 



10 



Men'shchikov et al.: Close binary and circumbinary torus of the Red Rectangle 



a radius of 430 AU (where dM/dr and M change their 
sign, Fig. ||), then the peak densities and mass-loss rates 
we derive would be lowered by a factor of 430 to the val- 
ues of riH ~ 10 10 cm -3 and M ~ 10~ 4 Mq yr _1 . It seems 
plausible that the hydrodynamical interplay between the 
gravitation, gas and radiation pressure, and bipolar out- 
flow have led to the present-day density distribution in- 
side the circumbinary torus and within the bipolar outflow 



cavities (see also Sect. 4.4). 



3.4. Near-infrared images 

In Fig. |^, our model images are compared side by side with 
our high-resolution images of the Red Rectangle in the H, 



K, Ks, and L c bands (Paper II; Paper III). The model in- 
tensity distributions at these wavelengths were convolved 
with circular Gaussian point-spread functions (PSF) of 
half-maximum widths of 47-76 mas corresponding to the 
resolutions of the reconstructed speckle images. In Fig. fj], 
we also present the model images at 0.61 fjxa and 10/xm, 
which can be qualitatively compared to an optical RG 610 
photog raph (|Cohcn et al. 1975| ), a deep coronagra phic H a 



image (Van Winckcl 2001 ), and an iV-band image (Waters 



ct al. 199S| ). The effects of saturation from long exposures 
in the optical images was roughly simulated by convolv- 
ing the model intensity distribution with a circular PSF 
of 8", whereas the 10 /xm model image was convolved with 
a PSF of 0'/7 which takes into account probable seeing 
conditions during the observations. Missing details on the 
PSFs and intensity scales in the papers mentioned above 
prevented us from using the images as constraints in the 
modeling and from making more accurate, quantitative 
comparisons. 

Two well-resolved bright lobes above and below the in- 
clined torus of the Red Rectangle arc seen in the observed 
images of Fig. [| along with the X-shaped spikes originat- 
ing deep inside the biconical outflow cavities. The spikes 
are bright enough to contribute to the intensity distribu- 
tion of the bright lobes (outflow cavities), making them 
appear broadened and even double-peaked in the highest- 
resolution Ks image (see also Sect. |3.6| ). The X-shaped 
spikes are caused most likely by a combination of the limb 
brightening and scattering by dust in the dense bound- 
aries of the biconical outflow cavities, as first indicated by 
polarization maps ( Perkins et al. 1981 ). Unidentified opti- 
cal emission bands, extended red emission, and the 3.3 /im 
PAH emission line also or iginate close to th e biconical sur- 
faces along the X spikes (Kerr et al. 1999). 



Extending our previous modeling of Paper II, in this 
work we attempted to simulate the X spikes appearing in 
the conical walls of the outflow cavities along the planes 
tangential to their surfaces. To accurately and consistently 
model the spikes, one needs a more sophisticated radia- 
tive transfer code than we have; moreover, the properties 
of dust and other species producing the X shape are very 
poorly known. For these reasons, we did not try to in- 
corporate the spikes in our 2D code in a self-consistent 



way as part of the radiative transfer iterations; instead, 
we decided to simulate the spikes for an illustrative pur- 
pose when computing the observable model fluxes. The 
spikes in the model images were produced by artificially 
increasing the scattering term of the source function in the 
cavities by 20 % between the plane tangential to the cavity 
surface and another plane closer by 5° to the symmetry 
axis. 

The spikes are clearly seen in all images of Fig. §, both 
observed and modeled, traced along the surface of the bi- 
conical outflow cavities to within the distances of 0'.'4 from 
the center. It is clear that the length of the spikes de- 
pends on the density distribution and dust properties, as 
well as on the sensitivity and resolution of observations. 
A higher signal-to-noise composite near-IR image of the 
Red Rectangle (Paper III) traces them to within 0'.'8 from 
the central binary. The spikes can be approximated by 
straight lines passing through the east and west peaks of 
the split northern and southern lobes and intersecting at 
the position of the central binary. The opening angle of 
the conical outflow cavities, measured between the spikes, 
is to «50°. 

The near-IR images of the Red Rectangle (Fig. 
clearly show the bright lobes which our model associates 
with a geometrically and optically thick torus inclined 
toward us in the south by ~11°. The images exhibit a 
compact, highly symmetric brightness distribution, very 
similar in shape and extent between the H and Lpah 
bands. As is becoming more apparent with highest res- 
olution of the K§ image, the bright lobes are divided by 
a very dark narrow lane, having a similar width in all the 
images. Even at the longest wavelengths there are no signs 
of either the direct light from the central binary, or an in- 
creasing contribution of hot emission (T<j ~ 1000 K) from 
the inner dust boundary of the torus. As discussed also in 



Paper III, such an appearance suggests a strong concen- 



tration of circumstellar matter toward the central binary 
in an axially-symmetric density distribution and a high 
gray optical depth due to large dust particles (Sect. fT6| ). 

3.5. Near-infrared color images 

High-resolution color images provide a sensitive measure 
of the relative intensity distribution at selected wave- 
lengths and a comparison of the observed and model color 
images gives a better idea of how closely the model resem- 
bles reality. Such a comparison is given in Fig. [?], where 
we display side by side the H~K, H—L c , and Ks~L c 
color images of the Red Rectangle with the corresponding 
model images at the same wavelengths. The correspond- 
ing near-IR speckle images (with resolutions of 68-7 6 mas) 
used t o produce the color images were presented in Paper 
H and Paper HI . 

The images in Fig. |?] exhibit a broad, relatively flat red 
plateau roughly 200 mas in size, covering the entire region 
of the bright lobes and the dark lane between them. The 
model color distributions in the near-IR bands are clearly 
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Fig. 6. Comparison of the H, K (76 mas resolution, Paper II ) and Ks, L c (47 and 64mas resolutions, Paper III ) 
speckle images of the Red Rectangle (upper row) with our model images at the same wavelengths (lower row). The 
model images were convolved with Gaussian PSFs having full widths at half-maximum equal to the resolutions of the 
corresponding speckle images. North is up and east is to the left in all panels. 




Fig. 7. Comparison of the H—K (76 mas resolution), H—L c (76 mas), and K$—L c (68 mas) color images of the Red 
Rectangle (3 upper left panels) with the corresponding model images (3 lower left panels). The color images were 
computed from only those parts of the respective images, where intensities are higher than 5 % of the peak. Mode l 
image of the Red Rect angle at 0.61 /iin (up per right panel) can be compared to a photograph in Cohen ct al. (1975 ) 
and to an H Q image in Van Winckel (2001). Model ima ge at 10 /im (lower right panel) convolved with a PSF of 0'.'7 
can be compared to an N image in Waters et al. (1998 ). North is up and east is to the left in all panels. 
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consistent with observations. Quantitatively, the reddest 
model colors of 3.0 mag, 5.4 mag, and 2.5 mag (in H—K, 
H—L c , and Ko,—L c ) are similar to those of 2.0 mag, 4.5 
mag, and 2.8 mag in the observed images. As the color 
images were computed from only those parts of the re- 
spective H, K, Kg, and L c images, where intensities are 
higher than 5 % of the peaks, the shapes of the outer edges 
of the color distributions essentially outline the brightest 
parts of the images. The blue X-shaped spikes are also 



clearly delineated in the color images. The red plateau 
marks the position of the very dense, optically thick, in- 
clined torus. The bluer colors outside the dense torus indi- 
cate optically thin regions in and around the outflow cav- 
ities which scatter the hot radiation of the central stars 
toward an observer. 



3.6. Intensity profiles 

Figure || presents a quantitative comparison of the model 
images, with available high-resolution optical and near-IR 



though the density distribution predicted by our model is 
more complex than a simple p oc r~ 2 profile assumed in 
their calculations. The two models have derived a similar 
bipolar geometry of the dusty envelope with density dis- 
tributions (almost) independent of the polar angle up to a 
limiting latitude, where the density drops to zero (or very 
small values) inside the biconical cavities. The authors 
assumed a full opening angle of 70° for the cavities, de- 
rived from low-resolution optical images (e.g., |Cohen et al 



1975). In our model, a smaller opening angle u> = 50° was 



adopted, which is consistent with our new high-resolution 
images. 



The model constructed by Lopez et al. (1997) differs 
from our model, in part due to the factor of 4 greater res- 
olutions of our new speckle images containing much more 
detailed information about the physical properties of the 
nebula. The differences can also be attributed to more 
extensive exploration of the parameter space in our mod- 
eling. The overall geometry of their envelope with bipolar 
cavities is very similar to the one adopted in Yusef-Zadeh 



imagesl of the Red Rectangle in terms of the intcnshy pro- 
files in two directions, parallel (PA = 12°) and perpendic- 
ular ( RA — 102°) to the projected symmetry axis of the 



et al. (1984) and in the present modeling, although the 



density is decreasing almost linearly with the latitude. The 

" 2 adopted by Lopez et al 



nebulal T he profi les (intensit}r cuts) were tak en from the 
images o fjPaper l| ( RG656, RG780), |Paper Il| (RG 715, H, 
K), and |Paper III] (Kg, L c , L PAH ). 

An inspection of the intensity profiles shows that the 
model is in good agreement with the observed intensity 
distributions of the Red Rectangle. There are almost no 
changes in the intensity profiles over all the near-IR wave- 
lengths except for those caused by slightly different resolu- 
tions. Such an invariance of the intensity profiles indicates 
that very large grains dominate the scattered and emitted 
radiation in the very dense circumbinary torus of the Red 
Rectangle. 

There are considerable discrepancies between the mo- 
del and observed images at optical wavelengths, where the 
bright lobes start to appear more diffuse and extended 
in the direction perpendicular to the symmetry (outflow) 
axis (see also images in Paper I ; Paper II ). In part, these 
deviations may be due to the approximations we used to 
describe the model geometry. The model density distribu- 
tion is assumed to be independent of the polar angle and 
this assumption might become more important at short 
wavelengths. We believe, however, that this morphological 
change can be naturally explained by additional contribu- 
tion of the extended red emission, known to be present 
in the Red Rectangle (but not in our model) at the op- 
tical wavelengths (0.55-0.75 /im). This emission veils the 
X-shaped spikes, making them appear relatively fainter in 
the visible than at the near-IR wavelengths, where the 
extended red emission does not affect the images. 

3. 7. Comparison with previous models 

Our results are consistent with the optical images of the 



radial density distribution p oc r 
(1997|) (with an order-of-magnitude jump upward beyond 



~ 100 AU) is relatively flat, containing most of the mass 
at the outer boundary. This would have a pronounced ef- 
fect (by a factor of ~ 2-3) on the beam-matched fluxes 
at long wavelengths, which their fit of the observed SED 
did not take into account. Very low densities of the inner- 
most regions in the model are inconsistent with the bipo- 
lar appearance of the Red Rectangle in the high-resolution 
near-IR speckle images shown in Fig. |^. 

Our previous model (Paper II) is very similar to the 
new model presented in this paper. As explained in Sect. [I], 
the old model was unable to reproduce the shape of the 
bright lobes of the Red Rectangle in the new Keck tele- 
scope images at the wavelengths of 3.1-3.3 /im. The model 
was also not very successful in fitting the observed SED 
over the UV to optical wavelength range and in explain- 
ing the flat distribution of the H—K color image over the 
bright lobes. All this stimulated us to reconsider all ex- 
isting observations and construct a more realistic model 
described in this paper. 

4. Discussion 

Below we interpret the results of our modeling of the Red 
Rectangle and discuss the properties and evolution of its 
close binary and dense circumbinary environment. 

4.1. Physical parameters of the binary 

The radio emission of a small H II region observed in the 
Red Rectangle (cf. Sects. 2.3, 3.2) indicates that there is a 
hot source of ionizing photons in the nebula. However, IUE 



Red Rectangle simulated by Yusef-Zadeh et al. (1984), al- 



spectroscopy (Sitko 1983) shows no signature of any hot 
T ii 2.5 x 10 4 K companion in the far-UV spectrum of the 
Red Rectangle, where one would expect to see the peak 
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Fig. 8. Comparison of the model intensity profiles with those from our optical and near-IR speckle images of the Red 
Rectangle. Also shown are the profiles of the point-spread functions and true (unconvolved) model intensities. Vertical 
dotted lines enable easier comparison of the relative offsets of the bright peaks (75 mas in the visible, 72 mas in H and 
K, 64.5 mas in JsTg, 65 mas in L c , and 66.5 mas in Lpah) from the center of the nebula. 
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of its radiation. This fact suggests that the hot invisible 
component has much lower luminosity than L+ « 6050 Lq 
of the observed post-AGB star. Taking into account in- 
terstellar extinction in our model (see Sect. 3.2), we can 
estimate that the well-observed SED sets an upper limit 
for the hot companion's luminosity of ~ 100 Lq. 

To illustrate this, in Fig. || we have added a hot (6 x 10 4 
K) blackbody to the stellar continuum of the post-AGB 
object, normalized so as to maximize its luminosity and 
yet to preserve our model fit of the observed SED. It is 
this hot additional component that produces rising far- 
UV fluxes toward the shortest wavelengths (A ^ 0.14/xm) 
in the stellar continuum and in the model SED with no 
interstellar reddening. Although the increasing interstel- 
lar extinction at these wavelengths (3 mag at 0.12/im) 
forces the model fluxes to decline, a bump in the SED due 
to the hot companion remains (Fig. ||). Our normaliza- 
tion implies that the invisible hot component's luminosity 
is a factor of at least ~ 60 less luminous than the post- 
AGB object, giving the upper limit of ~ 100 Lq. Together 
with the inferred effective temperature of T ^ 2.5 x 10 4 K, 
this limit marks an area in the Hcrtzsprung-Russell (HR) 
diagram populat ed by cooling white dwarfs with masses 
M > 0.23 Mq flDriebe et al. 199S| ). This can be inter- 



Table 2. Orbital elements of the close binary HD 44179. 



preted as the evidence that the close binary indeed consists 
of a post-AGB object and a relatively hot white dwarf. 
This idea is further supported by the spectroscopi c or- 

bital ejements nf t/hp binary's nhsprvpd cnmponpnt. ^/aXL 



Winck( |lLet al (1 QQr»|) Hpt.prminpH thp nrhital p eriod P = 



298 da ys an d the eccentricity e = 0.4 5, whereas [Waelkens 
et al. ( |l996| ) gave P = 318 days, and |Waters et al. (1998|) 
cited e = 0.38. We have recomputed the full set o f the 
latest orbital elements from original measurements ( Var| 
Winck el et al. 1995 ; Waelkens et al. 1996 ) using the com- 
puter code of Tokovinin (1992| )" 

The new orbital elements of HD 44179 are presented in 
Table |. With the viewing angle 9 V w 11° (Fig. 0) and the 
current mass of the post-AGB star of M* w 0.57 Mq (see 
for details Sect. the mass function f(M) = 0.049 M 
implies that the invisible companion is a low-mass object 
with AiyvD ~ 0.35 Mq. If one adopts a lower-limit effec- 
tive temperature of Lwd ~ 2.5 x 10 4 K (Sect. 2.3), one 
gets from the evolutionary tracks of white dwarfs ( Driebe 
et al. I 998) ) a very low luminosity of Lwd ~ 0.3 Lq for 
the st ar. In this paper we assume, however, that the hot 



Parameter 
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Value 


Error 


Units 


Period 


P 


322 


2.0 


days 


Periastron (JD 244+) 


T 


8075 


7.0 


days 


System velocity 


V 


19.8 


0.4 


kms -1 


Semiamplitude 


Ki 


12.2 


0.5 


kms -1 


Eccentricity 


e 


0.37 


0.03 




Semimajor axis 


a± sini 


0.34 


0.01 


AU 


Mass function 


f(M) 


0.049 


0.006 





the radiation transfer effects, with optical depths along 
different lines of sight greatly varying in such an axisym- 
mctric density distribution. The narrow dark lane appears 
quite naturally in bipolar envelopes like the one shown in 
Fig. [I] even when the "torus" density is independent of 
the polar angle. This is demonstrated by the model im- 
ages in Figs. ||, which show very good agreement with 
observations. 

A thin flat disk observed edge-on would appear as an 
ellipsoidal structure highly elongated at PA ss 102° (by a 
factor of ~ 5) perpendicular to the symmetry axis of the 
nebula, with a narrow and long dark lane (see, e.g., L-band 
model images of a circumstellar disk in suttner & Yorke 
2001). The expected appearance of a thin disk sharply 



contrasts the observed bipolar brightness distribution ex- 
tended along the symmetry axis. This rules out the flat 
disk geometry as a model for the interpretation of the ob- 
served properties of the Red Rectangle and other similar 
bipolar nebulae. Plausible distributions of dusty material 
may include either a strongly flaring, geometrically very 
thick disk or a toroidal structure (envelope) with biconical 
outflow cavities. 

We believe that it would be more accurate to describe 
the circumstellar environment of the Red Rectangle as an 
envelope with conical cavities, since the term "disk" is 
usually associated with a thin flat geometry. Many such 
dense envelopes are known among both young and evolved 
low-mass stars (L1551 I RS 5, HL Tau, IRC +10 216 are 
just a few of them; e.g., Men'shchikov fc Henning 1997 ; 
Men'shchikov et al. 19991 ; |White et al. 2000| ; |Men'shchikov 



compo nent's luminosity is iust at the upper limit of Lwn 
w 100 Lq and, therefore, significantly hotter, with Iwd 
« 6 x 10 4 K. 



4.2. The bipolar geometry 

Despite the narrowness of the dark lane, the circumbinary 
material we see in the images is definitely not confined to a 



flat disk. As discussed by Men'shchikov & Henning (2000) 



and Men'shchikov et al. (2001), the appearance of opti- 
cally thick structures with bipolar outflow cavities may 
easily be misleading. The narrow dark lane is a result of 



et al. 2001). In young stellar objects this geometry is natu- 
rally produced by the (radial) collapse of an initially spher- 
ical protostar, which eventually creates powerful jets and 
excavates bipolar (conical) outflows. In old stars this ge- 
ometry is naturally produced by the high (radial) mass 
loss they experience on the top of the AGB and by the 
fast bipolar winds expected at the later stages of their 
evolution. Close binarity of the stars plays an important 
role in the creation of such extremely dense, opaque bipo- 
lar envelopes, and the Red Rectangle is an outstanding 
example of such a structure. 
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4.3. The gas component 



4.4. Stability of the thick torus 



Although our model suggests a self-consistent picture of 
the Red Rectangle explaining a variety of observational 
data, the modeling was done using only the constraints 
from the dust continuum radiation (Sect. 2.2). Since the 
modeling directly predicts only the properties of the bi- 
nary and its bipolar dusty envelope, discussion of a wealth 
of information on the molecular and atomic species ob- 
tained during the past decades of observations is beyond 
the scope of this paper. Below we briefly discuss, however, 
a few aspects related to the circumbinary gas component. 

Where is the H II region located in the Red Rectangle? 
Our model requires such high densities in central regions 
of the nebula (n H ~ 2.5 x 10 12 cm- 3 , Sect. |J) that the 
relatively low-luminosity white dwarf would not be able 
to ionize hydrogen within a radius of ~ 17-170 AU, as 



seems to be required by the observed radio fluxes (Jura 
et al. ]L997 ) . This is true, however, only for a spherical 
distribution of circumstellar matter. The Red Rectangle 
has biconical outflow cavities with much lower densities 
of - 5 x 10 5 cm" 3 (Fig. |l|, S ect. fD|). A n estimate of the 
Stromgrcn sphere radius sq (Lang 1974) for the parame- 



ters of the hot white dwarf (Sect. 4.1) and electron tem- 
per at urecn^lFKTgrvesst^ 

the H II region fills in the outflow cavity of the dense cir- 
cumbinary torus of the Red Rectangle. 



Op tical spectroscopy of the Red Rectangle by Jura 
et al. (1997) revealed a two-component Ha emission con- 



sisting of a narrow spike of a full width of 20 kms at 
half-maximum, on top of a broad plateau with a width of 
200 kms -1 at zero intensity. Their analysis attributed the 
slow gas to the extended HII region, whereas the fast- 
moving ionized gas is likely to be confined to a much 
smaller zone that fits into the the binary's orbit. The 
ultraviolet spectra of CO and CI presented by plinski 



et al. ( L997) displayed narrow emission lines and wide ab- 
sorption bands with full widths of 25-30 kms -1 and 230 
kms -1 , respectively. The narrow emission and wide ab- 
sorption lines of CI seem to be related to the center of 
mass of the binary, although the absorption could be as- 
signed to a binary component, too. They noted a similarity 
of the lines with the two-component Ha profile observed 
by Jura et al. (1997 ), suggesting that the respective nar- 



row and wide lines arise in the same distinct environments. 

Clearly, new observations are needed to verify these 
findings of the high- velocity gas in the Red Rectangle and 
to place more precise constraints on the spatial distribu- 
tion of the observed species, before any realistic model 
could be able to interpret all existing data. Without any 
possibility to explicitly and self-consistently incorporate 
the gas components in our dust continuum radiative trans- 
fer modeling, there is no way to extract more information 
on the gases in the Red Rectangle from the model pre- 
sented in this paper. 



What keeps the gas and dust components of our model 
at high latitudes above the midplane, if the configuration 
is considered as stationary and long-lived, as proposed by 
[Jura et al. (1995| )? The torus-like (Fig. |l|) dense circumbi- 
nary envelope (Fig. |]) of the Red Rectangle does not re- 
semble a flat disk (Sect. [T^) and it may not necessar- 
ily be a stationary, centrifugally-supported configuration. 
Our model implies that it is a massive, selfgravitating, 
slowly evolving envelope ejected by the binary compo- 
nents, in which cavities were carved out presumably by 
a faster bipolar outflow. 

The idea of a stable disk in the Red Rectangle was in- 
voked by Waters et al. (199*2] ) to explain an unusual chem- 
ical composition of the spectroscopically observed post- 



AGB star (Van Winckel et al. 1995) that suggested sep- 
aration of the chemical elements due to the formation of 
dust and reaccretion of the depleted gas back onto the 



photosphere (Mathis & Lamers 1992). There is no good 
reason, however, why the circumbinary material should be 
necessarily distributed in a thin disk for this mechanism, 
which is based just on the radiation pressure on dust, to 



work. A long-lived disk was also proposed by Jura et al 



(1995 ) to explain unusual properties of dust in the Red 



Rectangle and possibly growth of dust particles to very 
large sizes of at least 200 /im. In a (stationary) disk, the 
dust grains would have much more time to grow than in 
the transient conditions of a normal wind of an AGB star, 
whose density p oc r~ 2 quickly drops as the envelope ex- 
pands. 

The growth rates due to coagulation of grains and ac- 
cretion of gas atoms on the grain surface is proportional 
to the rate of collisions and, therefore, to the gas and dust 
densities. As an example, the time scale (in years) for a 
grain to grow by accretion to a radius a can be written as 

* r><»( 77-r— ]!— If _](_H^(2) 



0.2 cm 



T 



\2 g cm 3 / \ rtH 

where p gr is the grain material density, jih is the number 
density of H atoms, T is the kinetic gas temperature, and 
a sticking probability of 0.1 and an abundance of 10 -3 of 
the accreted atoms were adopted. The high densities in 
the torus of the Red Rectangle (Sect. [O]) imply i acc ~ 60 
years for large dust particles to grow by accretion. This 
illustrates the possibility of a very rapid coagulation and 
growth of grains, on a crossing time scale of the dense 
toroidal envelope (~ 100 years, Fig. ||). Thus, there seems 
to be no need in a long-lived, rotating disk in the Red 
Rectangle to grow millimeter-sized dust particles. 

Although the dense environment is probably not sup- 
ported by the centrifugal forces, the toroidal envelope may 
well be quasi-static, balanced against the gravity by the 
gas pressure and by the radiation pressure on dust grains. 
The free-fall time scale for the gas component can be es- 
timated from 



t s « 118 



V100AU 



3/2 



2.12 M, 



1/2 



yr, 



(3) 
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where R c is the outer radius of the dense torus and M c « 
M* + Mwd + M is the mass inside it (as before, M is the 
mass of the torus). Substituting the model values, we find 
a free-fall time tg ~ 120 years. The dynamical expansion 
time scale of the torus due to the gas pressure can be 
expressed as 

1/2 



cxp 



170 



( R c 



ViooAuyv T 



yr. 



(4) 



As before, assuming that the gas temperature is equal to 
the dust temperature of ~ 10 3 K, we find i cxp rj 170 years. 
The order-of-magnitude comparison suggests that the gas 
component may be in a quasi-hydrostatic equilibrium. 

Estimating the maximum radius of dust particles a max 
that can be removed by the radiation pressure of the cen- 
tral star, we get an upper limit of 

o^^H^Y^V, (5) 



\6050L Q j \M r j \ p g 

where M r is the total mass within the radius r. For the 
high luminosity of the Red Rectangle derived in this mod- 
eling, the condition for dust grains to be gravitationally 
bound requires particles with a i; 0.2 cm at the inner 
dust boundary (r = 14 AU), where M r « M* + Jl%D. 
Since M C /(M+ + Mwd) ~ 2.3, the same condition of 
gravitationally bound grains at r — R c requires particles 
with a <; 0.09 cm. In our model, the radiation pressure 
cannot remove dust particles with a — 0.2 cm from the 
dense torus. The gravitation and radiation pressure on 
dust grains seem to balance each other in the envelope, 
causing its very slow quasi-static evolution. 

The dust grains outside the dense torus (r > 100 AU) 
have a normal size distribution in the model, therefore 
they are more efficiently accelerated by the radiation pres- 
sure and can drive a bipolar outflow. This can explain the 
observation by Jura et al. (1995 ) that the CO emission 
in the Red Rectangle is associated with the center of the 
mass of the binary and that it resembles a bipolar outflow 
with a speed of ~ 7 kms -1 . 

4.5. Evolution of the binary 

The presence of a large amount of matter in the close 
vicinity of the binary nucleus of the Red Rectangle sug- 
gests that the origin of the nebula may be associated with 
an unstable mass loss which accompanied formation of the 
post-AGB star, as noticed by Osterbart et al. (1997 ). This 
"unstable mass loss" may be related to the ejection of a 
common envelope in a close binary (Bond & Livio 1990; 
Yungelson et al. 1993). Thus we may consider an evolu- 
tionary scenario in which one component of a binary first 
forms a low-mass He or CO white dwarf, while another 
becomes an AGB star and then a CO dwarf, ejecting the 
currently observed envelope. In our analysis we followed 



the general outlines of the work of Nclcmans et al. (2000) 
on reconstruction of the evolution of close binary white 
dwarfs. For the evaluation of the lifetimes, radii, luminosi- 
ties, and core masses of the AGB stars we used the SSE 



package (Hurley et al. 2000). These quantities slightly de- 
pend on the specific prescriptions for mass loss adopted 
by [Hurley et al. (20001 ). 

In our discussion, the following notation is adopted. 
The initial masses of the components (on the main se- 
quence) are Mio and M 2 o (Mio > M 2 o) and their initial 
separation is ag. The mass of the primary after Roche lobe 
overflow (RLOF) is Mwdi and the separation of the com- 
ponents is a%. The mass of the secondary prior to RLOF is 
M 2 , its radius is R2, and the components' separation prior 
to the overflow is a 2 . After the RLOF, the mass of the 
secondary component is A%D2 (we consider the present 
post-AGB star as a nascent white dwarf) and the final 
separation of the components is af . We attempted to re- 
construct the formation of the binary of the Red Rectangle 
starting from the present configuration and "evolving" the 
system back in time toward the main sequence, to derive 
the initial parameters of the components. 

Changes of the components' separation in the com- 
mon envelope (CE) formed after the RLOF by the AGB 
star may be described by the equation of the balance of 
the binding energy of the giant's envelope and the differ- 
ence of the orbital energies of the pre- and post-CE binary 
( |Webbink 19841 ): 

M 2 (M^-M WD2 ) f M WD1 M WD2 M WD1 M 2 \ 
XR2 = a { 2^ 2aV~ J' (6) 

Here A is a dimensionless factor which allows one to ex- 
press the binding energy of the envelope as a function of 
Mg, M WD i, and R 2 ( |de Kool 1990| ). It depends on the 
structure of the stellar envelope, being a function of the 



evolutionary stage of the donor (see Dewi & Tauris 2000 
[Tauris fc Dewi 2001 ). The parameter a represents the ef- 



ficiency of the orbital energy deposition into the common 
envelope, depending also on the structure of the donor and 
the nature of the accretor. At the present state of the art 
in stellar evolution theory, one cannot disentangle these 
two parameters and it is reasonable to treat the product 
aX as a single parameter. 

The effective temperature T* w 7800 K of the post- 
AGB object suggests that it is now on the horizontal part 
of the white dwarf cooling track, therefore we can assume 
that the present L+ ks 6000 Lq is close to its luminosity 
before the CE ejection. Using "tracks" constructed by the 
SSE package, we estimate the mass M 2 , core mass M c , and 
radius R2 of the star at the moment when its luminosity 
is equal to the present value. We make then the usual 
assumption that Mwd2 = M c . The analysis shows that 
the initial mass M 2 o of the main-sequence progenitor of 
the post-AGB object had to be greater than 1.5 Mq, since 
stars with lower masses do not reach L± — 6000 Lq . 

One can also obtain an upper limit for the initial mass 
of the post-AGB object. Since the r.h.s. of Eq. (^) must 
be positive, we obtain the condition a 2 > (Mj/Mwra) a f- 
Assuming that the present mean separation of the com- 
ponents in the Red Rectangle represents the post-CE 
value of aj, we derive M 2 o 55 1.9 Mq. The narrow range 
1.5 ^ M 20 /Mq Ss 1.9 of possible initial masses for the 
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1.5 1.6 1.7 1.8 1.9 

Initial secondary mass (M a ) 

Fig. 9. The range of possible initial separations an of the 
components in the progenitor of the Red Rectangle binary 
depending on the initial mass M20 of the secondary com- 
ponent (the progenitor of the present post-AGB star). For 
the given M 2 n, the initial separation depends on the pa- 
rameter 7 in Eq. (|J). The upper edge of the shaded region 
corresponds to the initial mass M10 = 2.3 Mq of the pri- 
mary component, whereas the lower edge corresponds to 
equal masses of the components. 



present post-AGB object implies that now its mass is 
A%D2 ~ 0.57 Mq. Then the mass function and the in- 
clination of the orbit provide the mass A/wdi ~ 0.35 Mq 
of the first-formed white dwarf. 

The secondary component lost an amount of mass 
M 2 o — Mj as stellar wind, predominantly at the AGB 
stage. This widened the orbit, as the product of the sepa- 
ration and total mass was conserved. At the onset of the 
RLOF, the stellar radius R2 was equal to the radius of 
the Roche lobe, for which one can write (e.g., Paczynski 
19670 



i? L ~ 0.46 a 



m d 



1/3 



(7) 



where a is the separation of components, and m a 
are the masses of the donor and accretor, respectively. 
Substituting i?2, Mj, and Mwdi in Eq. ([?]), one can ob- 
tain an estimate of the separation of the components a 2 « 
670-5 90 R G (3.1-2.7AU) prior to the RLOF, for the ini- 



tial ms sses of M 



20 



1.5-1.9 Mr. 



q_ 



The mass Mwdi ~ 0.35 Mq of the first white dwarf 
suggests that it is a helium white dwarf descending from 
a red giant with M10 2.3 Mq, which had a degener- 
ate helium core. The radii of such stars depend only on 



the core mass (Ftefsdal & Weigert 1970). As for the sec- 
ondary, we assume that the dwarf's mass Mwdi is equal 
to the mass of the helium core of its progenitor. Then one 



can estimate the radius of the primary at the onset of the 
RLOF using, e.g., the core mass - stellar radius relation 



for giants of solar chemical composition (Ibcn & Tutukov 

10 35 (M 
ant with a 0.35 Mq 



wdi/M q ) . The radius of the gi- 
helium core was R\ ~ 47.5 Rq and 



the star should have had a deep convective envelope (more 
than 50 % of the mass in the outer convective zone, see the 
models of [Mengel et al. 1979|). Then the RLOF results in 



an unstable mass loss on a time scale intermediate between 
the thermal and dynamical ones. Nelemans et al. (20()()| ) 
have shown that, for comparable masses of the donor and 
accretor, the unstable loss of the stellar envelope most 
probably happens without spiral-in. 

The effect of the mass loss can be described in terms 
of the balance of the angu lar momentum. Following 
Paczynski fc Ziolkowski (1967 ), one may assume that the 
losses of mass and orbital angular momentum are related 
in a linear way, Jq~J\ = ^yJoAM/Mto. Here Jo an d J\ 
are the orbital angular momenta of the binary before and 
after the mass loss, AM is the mass lost from the sys- 
tem, and M t o = M10 + M20 is the total initial mass of the 
binary. The separation of components changes as 



oi = / MwMw 
a Im W diM 2 



M 



WDl 



Ma 



M, 



to 



1-7 



AMY 



M, 



to 



(8) 



where M2 is the mass of the secondary after the first mass- 
loss episode. For the stars in the mass range under consid- 
eration, a difference of, e.g., only ~ 0.1 Mq in initial mass 
is sufficient for the primary to complete its evolution in 
the red giants branch, while the secondary is still a main- 
sequence star. Therefore we can expect that the amount of 
accreted matter in the first mass loss episode in the system 
was negligible, i.e. that AM k, M 10 — Mwdi, M 2 ~ M 2 n. 
Equation (||) allows us to find the separation of compo- 
nents an of the initial binary if 7 is known. On the other 
hand, an can be found for any M 2 n using the core mass - 
radius relation and Eq. (fj]), since M 2 n < M10 < 2.3 Mq. 
This enables us to estimate the range of the parameter 
7. The solutions were obtained for 7 = 1.0-1.3, which 
are marginally compat ible with the estimates derived by 
Nelemans et al. (200C ) for low-mass close double white 



dwarfs. The range of a for different M 2 n is shown in Fig. [j| 
Several cautionary notes should be made on the analy- 
sis. In our scenario, the combined parameter aX in Eq. ^ 
has to be ~ 10. This may suggest that the release of the 
internal thermodynamic energy may be involved in un- 
binding the envelope ( Han et al. 1994 , 1995 ). Although 
the values of the structural parameter A have never been 
derived for stars with M < 3 Mq , large values of A may be 



characteristic for the tip of the AGB (see Dewi & Tauri 
2000). In addition, on the AGB the effect of the orbital 



shrinkage may be somewhat offset by the wind outflow, as 
has been noticed by [ben (200C| ) . Strictly speaking, Eq. (JsJ) 
becomes then invalid, but effectively this would imply an 
increase of aA. 

Our estimate of Mwdi is compatible as well with an 
assumption that this dwarf is a descendant of a star with 
an initial mass of about 2.5-3 Mq, which overfilled the 
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Red Rectangle: Evolutionary scenario 

M 10 =2.3,M 20 =1.9, ao=126 
M 10 =2.3, M 2 =1.9, ao=126 
M WD1 =0.35, M 2 =1.9, a.,=546 



Main-sequence binary 
Red giant + MS star ^ <j 

Unstable loss of RG envelope — * 



He white dwarf + MS star 
He white dwarf + AGB star 
Outflow of common envelope 



M WD1 =0.35, M 2 =1.9, a.,=546 



o 



=0.35, M 2 =1 .74, 32=589 



M WD1 =0.35, M WD2 =0.57, a,=192 



Fig. 10. The sequence of evolutionary stages of the binary 
resulting in the formation of the Red Rectangle nebula. 
The binary parameters are given in solar units. 



Roche lobe prior to the helium ignition in the core and 
produced a helim star that later became a CO dwarf ( Ibcr| 
& Tutikov 1985|; |Han et al. 2000] ) . However, for the for- 



mation of a 0.35 Mq dwarf the RLOF had to occur when 
i?2 20 Rq and the resulting separation of components 
a2 would not be sufficient to harbor an AGB star with a 
0.57 Mq core. 

The RLOF occurs close to the tip of AGB and some 
mass at this evolutionary stage has already been lost by 
stellar wind. For the initial masses of 1.5-1.9 Mq of the 
post-AGB object in the Red Rectangle, this amounts to 
the loss of 0.35-0.16Mq. The mass of the common enve- 
lope is therefore constrained to M — 0.58-1.17Mq in our 
evolutionary scenario. This value derived from the above 
analysis allowed us to constrain the dust-to-gas mass ratio 
Pd/p in a way similar to the approach used in our model 
of IRC +10 216 ( [Men'shchikov ct al. 200l] ). In fact, from 
our continuum radiative transfer modeling we know the 
total mass of dust Md rs 0.015 Mq in the Red Rectangle, 
which is contained in the dense torus. The estimated range 
of possible masses of the circumbinary torus implies that 
Pd/p = 0.021-0.01, rather high values compared to the 
usual Pd/p ~ 0.003 in the winds of single mass-losing AGB 
stars. The higher value seems to be too high considering 
the abundances of elements that can make up dust grains. 
We therefore adopted in this study the lower end of the 
interval, Pd/ P = 0.01, which corresponds to the total torus 
mass M rs 1.2 Mq and the initial mass M20 ~ 1.9 Mq of 
the star. 

Figure ^ illustrates and summarizes the evolution- 
ary scenario presented above. The initial masses of the 
components in the progenitor binary system of the Red 



Rectangle are very similar to each other, but not equal. 
Whereas a 2.3 Mq star completes its evolution to the tip 
of the red-giant branch in about ~ 800 Myr, a 1.9 Mq star 
evolves to the tip of AGB in ~ 1600 Myr. Thus, there may 
be a difference of about 800 Myr in the ages of the binary 
components of the Red Rectangle. Within this time a he- 
lium white dwarf would cool down to ~ 0.01 Lq (Driebe 
ct al. 1998). If its present luminosity has indeed a much 
higher value of Lwd Si 100 Lq, this might mean that 
the He white dwarf has accreted some hydrogen from the 
AGB star wind and experienced one or several hydrogen- 
burning flashes. Then the latter brought it close to the 
high luminosity and temperature "knee" of the cooling 
track again. 

4.6. Eccentricity of the orbit 

It is well known that tidal interactions in close binaries 
tend to quickly circularize their orbits. The large eccentric- 
ity e = 0.37 of the close binary in the Red Rectangle has 



been interpreted in two different ways. Van Winckel et al 



(1995) suggested that increased mass loss, when the star 



passes periastron, may be the process responsible for the 



increased eccentricity, whereas in a later paper (Waelkens 
|ct al. 1996 ), the authors suggested that gravitational in- 
teraction of the binary with the circumbinary disk may 
be the reason. Recently, Sokcr (2000] ) argued that it is the 
mass loss at periastron which increases the binary eccen- 
tricity and that the inner disk mass is too small for the 
other mechanism to be efficient in this and other similar 
objects. 

In our analysis, the initial orbit of the binary is circu- 
lar. It is usually assumed that the unstable mass loss by 
Roche lobe filling stars does not produce eccentric orbits. 
However, the dynamical time scale for the AGB star was 
shorter than the orbital period of the binary and we may 
speculate that the relatively fast loss of the envelope made 
the orbit elliptic. Later, the eccentricity can increase as a 
result of the binary-disk interaction. 

Our modeling has shown that the torus' mass in the 
Red Rectangle may be large enough to account for the 
eccentricity changes. The mass in question resides inside 
a radius of r « 6 a\ and it must be sufficiently large for 
the circumbinary disk to efficiently drive the increasing 
eccentricity (Artymowicz et al. 1991). A lower bound on 
the model torus mass within this radius can be found 
by extrapolating the density at the inner dust boundary 
(4.2 x 10 -12 gem -3 ) into the dust-free zone. This gives a 
mass of ~ 0.001 Mq and a disk-to-binary mass ratio qd 
of the same order. Assuming an initial eccentricity of 0.1, 
wc find from e w 0.002 q d 2n/P QArtymowicz et al. 199l| ) 
a high rate of changes on a time scale of e/e ~ 7 x 10 3 yr. 
We note, however, that this estimate is just an illustration 
and that it remains far from being certain that this specific 
mechanism, assumptions, and formulae are applicable for 
the binary's environment existing in the Red Rectangle. 
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5. Conclusions 

Recent diffraction-limited near-IR speckle images of the 
Red Rectangle in the wavelength r ange 2.1- 3.3 /im with 
angular resolutions of 44-68 mas ( Paper III ) as well as 
the pre vious op t ical and near-IR speckle images at 0.7- 
2.2 /xm ( Paper I ; Paper II ) revealed a geometrically thick 
circumbinary torus with bipolar outflow cones and X- 
shaped spikes originating deep inside the cavities. This 
multiwavelength set of high-resolution images enabled us 
to reanalyze most of the existing observations using our 
two-dimensional radiative transfer code. Results of this 
study are summarized below. 

Distance. An important by-product of the modeling 
of the Red Rectangle is a new determination of its distance 
D « 710pc (with model uncertainties of about 10%), 
which is twice as large as the commonly adopted value 
of 330 pc (Sect. UM. The new distance, based on the ac- 



count for interstellar extinction in our model (Sects. 2.4, 
3.2), is consistent with high luminosities expected from 
the stellar evolution theory. Based on this distance, we 
reconstructed physical parameters of the binary and its 
circumbinary torus (Table [l]) . 

Close binary. In our model, the observed compo- 
nent of the spectroscopic binary is a luminous, low-mass 
post-AGB star with M* « 0.57 M , T+ « 7750 K, and L± 
« 6050 Lq . It is a product of the mass loss by the sec- 
ondary component in a close binary system. We identified 
the now invisible descendant of the primary component 
with a relatively hot white dwarf with A%d ~ 0.35 Mq, 
Twd ~ 6 x 10 4 K, and Lwd ~ 100 Lq. The presence and 
parameters of the compact degenerate star _m the close 
binary HD 44179 have been deduced (Sect. 4.1) from (1) 
the spectroscopic mass function f(M) = 0.049 Mq, (2) 
presence of a compact H II region ionized by a hot source 
of radiation, (3) upper limit of the contribution of the hot 
object to the continuum in the far UV. 

Circumbinary structure. The intensity distribu- 
tion of the high-resolution images is definitely inconsis- 
tent with the flat disk geometry frequently used to visual- 
ize bipolar nebulae (Sect. 4.2). A geometrically very thick 
density distribution of a compact, dense torus with bicon- 
ical outflow cavities (Fig. 0) is best suitable for reproduc- 
ing the observed images. The opening angle of the cavities 
is 50° and the observer's viewing angle is 11° below the 
midplane. 

Density distribution. Although the nebula extends 
to at least i?2 ~ 4 x 10 4 AU from the central binary, 
most of its mass of M w 1.2 Mq is contained in the 
extremely dense torus having an outer radius of 100 AU 



(Sect. 3.3). The inner dust boundary of the torus is lo- 
cated at a distance R\ « 14 AU from the center, where 
gas densities reach values of p w 4.2 x 10~ 12 gem -3 
(riH ~ 2.5 x 10 12 cm -3 ). The density of the outflow cavi- 
ties is many orders of magnitude lower than the torus den- 
sity in the region dominated by the dense torus (Fig. |J), 
whereas in the outer regions (r > 800 AU) the outflow 
cones are denser than the rest of the toroidal envelope. 



Dust properties. Our model of dust in the circumbi- 
nary torus of the Red Rectangle has two distinct compo- 
nents whose chemical composition is not constrained by 
observations. Amorphous carbon grains with radii a in the 
range of 0.005-600 /im and typical interstellar size distri- 
bution dn/da oc a -3 5 exist mainly in the outer regions of 
the toroidal envelope (r > 100 AU). Most of the dust mass 
is contained in very large (a = 0.2 cm) particles of the mas- 
sive, dense torus, which produce an almost gray optical 
depth of r » 47. Although there exist more dust compo- 
nents somewhere in the nebula, including crystalline sil- 
icates and PAHs, their mass must be small compared to 
the mass of the very large particles, which dominate the 
observed appearance of the Red Rectangle. 

Evolution of the binary. Based on our estimates 
of the masses of the Red Rectangle binary components, 
their separation and luminosities, we suggest an evolu- 
tionary scenario for the formation of the nebula, in which 
components had initial masses of about 2.3 and 1.9 Mq 
and a separation of ~ 130 Rq. The scenario associates the 
formation of the Red Rectangle nebula with the ejection 
of the common envelope upon the Roche lobe overflow by 
the present post-AGB star. 
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